Frontiers in Molecular Biosciences

Inadequate autophagy promotes smoking-induced cellular senescence involved in chronic obstructive pulmonary disease (COPD) pathogenesis. Transcription factor EB (TFEB) is a master regulator of the autophagy-lysosome axis. For the first time, we investigated the therapeutic potential of pemafibrate, a putative TFEB inducer. COPD lung epithelial cells showed reduced TFEB expression. Pemafibrate enhanced autophagy/mitophagy flux and restored lysosomal acidification observed during cigarette smoke (CS) extract exposure in human bronchial epithelial cells, resulting in reduced cellular senescence. TFEB knockdown demonstrated involvement of pemafibrate-induced TFEB in these effects. Pemafibrate induced TFEB expression, mitigated alveolar enlargement and airflow obstruction, and attenuated the CS-induced increase in static lung compliance in a long-term CS-exposed mouse model. It reduced the CS exposure-induced cellular senescence, possibly through autophagy/mitophagy, as suggested by bulk RNA sequencing of mouse lungs. A retrospective cohort study showed that patients given pemafibrate displayed attenuated FEV1.0 decline compared with those given bezafibrate or fenofibrate. In conclusion, pemafibrate is a promising therapeutic agent for COPD, potentially exerting its effects through the regulation of the TFEB-autophagy/mitophagy-lysosome axis.

Next-generation sequencing (NGS) for the detection of somatic variants has become the method of choice in a variety of molecular oncology fields and in the clinic. Its use ranges from sequencing entire tumor genomes and transcriptomes to targeted clinical diagnostic gene panels. The NYU Langone Genome PACT (Profiling of Actionable Cancer Targets, LG-PACT) assay is a qualitative in vitro diagnostic test that uses targeted next generation sequencing (NGS) of formalin-fixed paraffin-embedded (FFPE) tumor tissue matched with normal specimens from patients to detect gene alterations in a targeted panel covering 606 genes and the TERT promoter. Indications for testing are cancer (solid tumors and hematological malignancies) where a mutational profile from multiple genes would be informative for disease stratification, prognosis, or treatment options including targeted therapies and eligibility for clinical trials. The test is intended to provide information on somatic mutations including point mutations, small insertions/deletions (indels), and copy number aberrations for diagnostic and treatment decisions. LG-PACT is a United States Food and Drug Administration (FDA) cleared diagnostic test (510K: K202304). The clinical interpretation of sequencing data of molecular tumor markers from NGS encompasses automated variant calling tools with human interpretation. This final mostly manual review of data step is intensive, involving highly trained scientists, encompassing literature review, interpretation and clinical tier classification by pathologists, who then provide a complete molecular diagnostic report to the treating oncologists. We provide analysis of 1339 clinical genomic profiles from 31 different cancers and their subtypes, comprising of central nervous system (CNS) 792 (59%) cases (incl. meningioma, glioma and glioblastoma), with 267 (20%) cases predominantly of lung, pancreatic and colorectal and 280 of others (21%). Here, we present the technical challenges of validating an NGS oncological diagnostic targeted assay for clinical grade accuracy and sensitivity for patient care. We show how copy number alterations provide a more comprehensive description of the tumors genomic profile. We then outline the utility of targeted panel sequencing based on certified pathologist selection of reportable variants for our current patient cohort. Where analysis of variant detection has led to 49.4% (661/1339) of our clinical tumor samples containing mutations in known therapy targeted genes, 35.6% (477/1339) with mutation detected in other genes, and 15% (201/1339) cases being negative.

Summary: Prot2Surf is a software tool designed for the characterization and prediction of protein association to surfaces. In this application note, Prot2Surf was tested using catalytic domains of the lytic polysaccharide monooxygenases (LPMOs), interacting with native surfaces. The results show that the software can efficiently analyze key binding features, including protein-surface distances, distances between catalytically reactive atoms, and the orientation angle between surface chains and the protein. These features are essential for distinguishing productive binding poses in these protein-surface systems and for understanding interaction patterns that provide guidance on protein engineering. Prot2Surf performs these analyses within seconds to a few minutes, providing a fast and accessible framework to post-process and characterize protein-surface encounter complexes. Availability and implementation: Prot2Surf, which is written in Fortran90, is documented and freely available as open source on GitHub: https://github.com/TUNNELING-GROUP/Prot2Surf. In order to run Prot2Surf, users should also install the SDA software package which is freely available at https://www.h-its.org/downloads/sda7/.

Antimicrobial resistance is a public health challenge, driving the need for rapid, cost-effective diagnostic support tools. Artificial intelligence (AI) may enable prediction of susceptibility to untested antibiotics from known susceptibility results, but prospective clinical validation is required before routine use. We evaluated an AI-based decision support method, trained on invasive isolates from the European Surveillance System (TESSy), for prediction of antibiotic susceptibility in clinical Escherichia coli urine isolates. The evaluation included 99 E. coli isolates from urine samples with diversity in age, sex, and antibiotic susceptibility. Predictions were evaluated for 14 antibiotics using patient metadata and susceptibility results for 4-8 antibiotics as input. Prediction uncertainty was handled using conformal prediction, allowing abstention when confidence was insufficient. EUCAST disk diffusion test results were used as reference and genomic sequence data was used to explore mechanisms of the AI performance. Without conformal prediction, 84% of predictions were correct when susceptibility results of six antibiotics were used to predict susceptibility to eight additional antibiotics. Across all predictions generated using susceptibility results for six antibiotics as input, the major and very major error rates were 19% and 12%, respectively. Prediction errors varied between antibiotics and were associated with certain phenotypic and genotypic resistance patterns. Conformal prediction reduced errors but increased abstentions; at confidence levels of 90%, 95%, and 97.5%, the model abstained in 9.6%, 14%, and 22% of instances. The method showed promising performance, but its clinical use remains limited and may require diagnostic data beyond susceptibility test results and demographic variables.

Aptamers are structured nucleic acid ligands capable of high affinity, high specificity molecular recognition generated using variations of the SELEX (Systematic Evolution of Ligands by Exponential Enrichment) process. However, SELEX often produces sequences that enrich yet may lack binding efficacy. We propose a measure called the Ruggedness Composite Index (RCI) along with a method for computing it, that can be used to distinguish binding-competent ('active') aptamers from weak or non-binding ('inactive') aptamers. Given a set of aptamers, RCI incorporates information on their fragmentation (landscape partitioning), basin entropy (metastable state distribution), cumulative density irregularity (non-uniform occupancy), and structural energy correlation length (structure-energy coupling scale). We test whether secondary-structure folding energy landscape topology distinguishes active from inactive aptamers using a multiscale level set framework across six datasets. Active aptamers show lower RCI values and occupy smoother, funnel-like conformational spaces, while inactive aptamers show higher RCI values, reflecting fragmented, high-entropy landscapes. By contrast, classical thermodynamic features, such as minimum free energy, show limited discrimination between active and inactive aptamers. In all datasets, sequences that exhibit enrichment which is not monotonic but lack specificity exhibit elevated ruggedness, indicating landscape topology can predict non-specific enrichment. These results indicate that folding landscape organization can be used as a predictor of aptamer activity and establish RCI as a simple, mechanistically interpretable measure for improving candidate prioritization, especially in therapeutic aptamer discovery.

Rationale: Pulmonary hypertension is defined by pulmonary hemodynamics, but diagnostic and prognostic biomarkers remain limited. Nuclear magnetic resonance (NMR) spectroscopy provides detailed insights, particularly in the lipid metabolism. Objectives: To explore circulating NMR-derived metabolites and lipoprotein-related parameters for their association with pulmonary hemodynamics and to analyse their prognostic properties in pulmonary arterial hypertension (PAH). Methods: Retrospective analysis of a PAH cohort with complete diagnostic workup including right heart catheterization and baseline serum samples, from the prospective GRaz Pulmonary Hypertension-Metabolism (GRAPH-M) registry. Measurements: NMR-derived metabolites and lipoprotein-related parameters were analyzed for their association with clinically relevant parameters of PAH. We defined PHIHDL, a score derived from high-density lipoprotein (HDL) related measures based on their strong association with pulmonary hemodynamics, and evaluated its prognostic value. Results: We included 100 patients with PAH treated at the PH clinic of LKH University Hospital, Medical University of Graz, between 2011 and 2021. Age was 61{+/-}15 years, female/male ratio 2.5, BMI 26 {+/-}7 kg/m2, mPAP 41{+/-}16 mmHg, PAWP 8.8{+/-}3.2 mmHg, PVR 8.0{+/-}4.9 WU, and median survival was 8.0 years. During follow-up, 46 patients died. We identified a cluster of 12 HDL-related measures that showed significant inverse association to pulmonary hemodynamics and derived PHIHDL from the reversed scaled average of these particles. PHIHDL was associated with all-cause mortality after adjustment for age and sex (HR 2.96, 95% CI 1.52-5.70), independent of the clinical risk scores COMPERA 2.0 and REVEAL Lite2. Conclusion: PHIHDL, a pulmonary hemodynamics-based metabolomic score, provides independent prognostic information beyond established risk scores in PAH.

Chemoradiotherapy (CRT) is the standard-of-care therapy for many solid malignancies, yet predictive biomarkers of treatment response remain limited. We identified a germline single nucleotide polymorphism (SNP) in an intrinsically disordered region of the lysine demethylase KDM3C/JMJD1C (p.S464T) that is associated with CRT outcomes in locally advanced rectal cancers (LARC) and head and neck squamous cell carcinoma (LA-HNSCC). In silico modeling with AlphaFold predicted S464T substitution influenced interaction between phosphorylated KDM3C and RNF8 FHA domain. In cellular models, conversion of S464 to T464 increased sensitivity to DNA-damaging agents. S464T substitution impaired damage-induced MDC1-RAP80 signaling and downstream RAP80-BRCA1 colocalization. SNP carrying cells impaired DNA repair causing genotoxic stress that is associated with increased cGAS-cGAMP innate immune signaling and increased apoptosis. Population analyses with the SNP highlighted an increase incidence of UV-induced skin and other cancers, linking inherited variation in the chromatin regulatory gene KDM3C to genome instability, cancer risk, and therapeutic vulnerability.

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.

Nanodiscs and bicelles are widely used as membrane mimetics for structural studies of membrane-associated systems. Studies have reported that their magnetic alignment behavior and phase stability are highly sensitive to composition and temperature. In this study, we systematically investigate the effects of cholesterol on bicelle formation and magnetic alignment in DMPC + 0.2 glycyrrhizic acid (GA) systems using a combined 31P and 14N solid-state NMR experimental and simulation-based approach. Temperature dependent 31P NMR spectra reveal a clear transition from vesicle dominant to aligned bicelles/nanodsics phase, while 1N quadrupolar splitting and lineshape analysis provides quantitative insights into heterogeneous lipid bilayer populations, distinguishing large aligned nanodiscs (B(L)), small nanodiscs (B(S)), and isotropic/random components (B(R)). A strong correlation is observed between the 31P derived bicelle fraction and the 14N B(L) population, confirming that macroscopic alignment in the presence of an external magnetic field directly reflects the growth of large, well-ordered nanodiscs. Cholesterol is found to play a critical dual role by modulating membrane order and curvature. At low cholesterol concentration (0 to 5 mole percent), nanodiscs alignment occurs gradually with increasing temperature, while at higher cholesterol concentration (15 to 25 mole percent), the alignment is delayed and accompanied by broader spectral features, indicating structural heterogeneity. Notably, 10 mole percent cholesterol consistently provides the optimal balance, enabling efficient temperature dependent conversion to aligned bicelles while maintaining high B(L) populations (about 70-80 percentage) and minimal isotropic fractions. In contrast, higher cholesterol maintains significant B(S) and B(R) populations, even at elevated temperature. The 14N quadrupolar coupling (Cq is approximately 8.5 to 9.2 kHz for aligned nanodiscs) remains nearly invariant across compositions, showing that cholesterol does not change local headgroup dynamics but instead redistributes lipid populations. These findings establish a combined 31P and 14N solid -state NMR approach provides a valuable platform for quantitatively correlating membrane structure, dynamics, and alignment, offering practical guidelines for optimizing bicelle systems for high resolution solid-state NMR studies of membrane associated biomolecules.

Background: The condition of a patient with acute stroke often changes within hours of ICU admission. Prognostic work here targets fixed endpoints predicted from admission data, and trajectory phenotyping assigns one label per patient. We used longitudinal ICU data to identify interpretable dynamic clinical states, characterize transitions between them, and relate the current state to later events. Methods: Retrospective cohort study of 6368 adults with acute stroke in MIMIC IV v3.1. The first 72 h were divided into twelve 6-hour windows, and a hidden Markov model was fitted to 21 neurological, physiological and organ support variables. State number was chosen against criteria fixed before fitting: statistical fit, restart stability, state occupancy and clinical interpretability. Generalized estimating equations related the current state to new mechanical ventilation and vasopressor use within 12 h, and to ICU death within 72 h. Eleven sensitivity analyses assessed the robustness of the state solution. Results: Four states were selected: neurologically preserved-low support, neurological impairment low support, impairment renal dysfunction and impairment-respiratory support (63.3%, 7.8%, 11.8% and 17.1% of windows). Within 72 h, 40.3% of patients changed state at least once, and transitions ran in both directions rather than along a single severity gradient. States were identified without outcome data, yet ICU mortality by last state ranged from 2.9% to 43.9%. Adjusted for age, sex, subtype and Charlson index, the current state remained associated with organ-support escalation and death. State prevalence differed by at most 1.1 percentage points between training and test sets, and 10 of 11 sensitivity analyses gave a stable four-state solution (ARI 0.754 0.955). Conclusions: The early ICU course of acute stroke can be represented as movement among a small number of clinically interpretable states. The representation was reproducible in a held out set and across admission eras, but requires validation in an independent database before any clinical use.

Biomolecular condensates formed by intrinsically disordered proteins require molecular models that accurately describe proteins in both dilute solution and condensed phases. Explicit-solvent coarse-grained models offer an attractive balance between chemical resolution and computational efficiency. Yet, it remains unclear whether improving dilute-state properties is sufficient to obtain an accurate description of condensates. Here, we address this question by introducing minimal modifications to the Martini 3 force field that combine recent advances in bonded interactions with refined protein-water interactions and strengthened glycine self-interactions, while preserving the underlying chemical transferability of the model. The resulting model substantially improves the description of single-chain conformations across a diverse benchmark of disordered proteins. We then investigate phase separation of the well-characterized low-complexity domain of heterogeneous nuclear ribonucleoprotein A1 and its sequence variants. The model reproduces several key physicochemical properties of biomolecular condensates, including chain expansion in the dense phase, sequence-dependent intermolecular contacts, protein diffusion and its relation to single-chain dimensions, and hydration, while revealing quantitative limitations in condensate density, phase equilibria, and ion partitioning. Our results show that improving dilute-state behaviour translates into a better description of condensed-phase properties, including condensate density, but is not sufficient to quantitatively reproduce the equilibrium between the dilute and dense phases.

Background: Higher Body Mass Index (BMI) is an established risk factor of sleep disturbance. It is not known if the effect is homogeneous across the lifecourse or if there is a particular time point in life that might be best to target. Methods: Two-sample Mendelian randomization (MR) was used to investigated the effect of childhood adiposity (adjusting on adulthood adiposity and obstructive sleep apnea (OSA)) on insomnia, morning chronotype, sleep duration, daytime sleepiness and daytime napping. Similarly, total, and direct effect of adulthood adiposity on these outcomes was explored. We used summary statistics from a genome-wide association study (GWAS) of UK Biobank for childhood and adulthood adiposity (n=453,169) and large-scale consortia of OSA (Million Veteran Program) (n=410,268), insomnia, and chronotype (23andMe) (n=1,978,022 and n=248,1000, respectively). Results: Two-sample univariable MR analysis provided no evidence of an effect of genetically predicted childhood adiposity on later life insomnia (Odds ratio (OR)= 0.94, 95% Confidence interval (CI)= 0.87, 1.03). Whereas, multivariable MR (adjusted for adulthood adiposity) analysis provide strong evidence of direct protective effect of genetically predicted childhood adiposity on later life insomnia (OR= 0.70, CI= 0.64, 0.77). Further, both in univariable and multivariable MR, a strong positive effect of increased childhood body size on morning chronotype was observed (OR= 1.16, CI= 1.01, 1.33 and OR= 1.36, CI= 1.15, 1.62, respectively) after accounting for adulthood body size. In both analysis the estimate did not change considerably after aditionally adjusting for OSA. However, childhood and adulthood adiposity found to be associated with OSA and OSA with insomnia. In both univariable and multivariable analysis, increased body size in adulthood increased the risk of having insomnia and a morning chronotype. Conclusions: The findings suggest that higher body size in childhood is not a risk factor for later life insomnia, whereas higher body size in adulthood was. Further, if healthy body size is maintained in adulthood, high childhood adiposity may decrease the risk of insomnia and increase the risk of being a morning person in later life. Keywords: childhood, adulthood, obesity, insomnia, morning chronotype, medelian randomization

Background: Oral dysbiosis has been associated with oral squamous cell carcinoma (OSCC); however, most microbiome studies rely on 16S ribosomal RNA (rRNA) gene sequencing, limiting species-level taxonomic resolution. Methods: Dental plaque, saliva, and tumor tissue samples from 10 patients with OSCC and dental plaque and saliva samples from 10 healthy controls were analyzed in this exploratory cross-sectional study. DNA was extracted and subjected to shotgun metagenomic sequencing using the Illumina MiSeq platform. Sequence reads were quality filtered with fastp, taxonomically classified using Kraken2 v2.1.3, and species-level abundances were re-estimated with Bracken v2.9 following the removal of human reads and low abundance taxa. Relative abundances were compared using the Mann Whitney U test with the Benjamini Hochberg false discovery rate correction, while the Bray Curtis principal coordinate analysis was used as an exploratory approach to visualize microbial community patterns. Results: Shotgun metagenomic sequencing revealed distinct bacterial community profiles across the oral microenvironment. Dental plaque exhibited the highest taxonomic diversity and relative abundance. The control plaque was enriched in Streptococcus koreensis, Capnocytophaga sp. oral taxon 878, Treponema sp. Marseille Q4132, and Leptotrichia sp. oral taxon 498, whereas the plaque from patients with OSCC showed a higher relative abundance of Pyramidobacter piscolens, Parvimonas parva, and Gemella sanguinis. Salivary samples displayed lower diversity and a more homogeneous composition, predominantly comprising Capnocytophaga endodontalis, Prevotella jejuni, Aggregatibacter aphrophilus, and Gemella sanguinis. The tumor tissue showed relatively higher abundance of Sellimonas catena, Escherichia coli, Solobacterium moorei, and Lacrimispora sp. HJ 01. Conclusions: This exploratory study provides species-level characterization of the oral microbiome across multiple oral microenvironments in OSCC and generates hypotheses for future integrative metagenomic and functional studies investigating the potential contribution of oral bacterial communities to OSCC pathogenesis.

Obesity is a significant public health concern. Early-onset obesity in the context of rare disease can reflect genetically-mediated pathology or elevated susceptibility through indirect mechanisms. Mapping the diverse characteristics and needs of young people with obesity in the rare disease population is a first step toward mechanistic and translational research. We carried out a retrospective comparative analysis of demographic, genotypic, phenotypic and health service utilisation data for young people with obesity (cases: n=500) and without obesity (controls: n=11,444) from the UK 100,000 Genomes Project rare disease cohort. Cases and controls were recruited prior to genomic diagnosis, across clinical disorder categories. We observed significant association between socioeconomic deprivation and obesity risk. Young people with obesity had significantly higher utilisations of acute care and mental health services, indicating an overall higher health burden. A curated panel of 519 candidate obesity-associated genes demonstrated aggregate association with obesity, although no single gene reached significance. Phenotypic comparison between cases and controls highlighted increased multi-organ and neurological system involvement, highlighting the overlap between neurodevelopmental and obesity risks. Within the case group, we conducted cluster analysis to identify early-onset obesity groups with different phenotypic profiles, potentially arising from different causal pathways - this identified six obesity subgroups of interest, with differing involvement of neurodevelopmental and other systems. Our study confirms that obesity co-occurs with a wide range of factors within the rare disease population, and is associated with significant physical and mental health needs, requiring holistic lifelong care.

Background: The Boston Criteria v2.0 represent the gold standard for diagnosing Cerebral Amyloid Angiopathy (CAA), but their application is currently precluded in mixed small vessel disease (SVD), where deep and lobar hemorrhages coexist. The aims of this study are: (i) to determine which cerebrospinal fluid (CSF) biomarker (A{beta}42, A{beta}40, A{beta}42/40 ratio) is the best candidate to support the CAA diagnosis; (ii) to define a data-driven cut-off, and (iii) to explore if a biomarker-integrated classification significantly improves the phenotypical concordance with the suspected predominant SVD (CAA vs. arteriosclerosis). Methods: We analyzed data from a retrospective multicenter cohort of patients with suspected CAA, defined as probable CAA (Boston criteria v2.0) but allowing deep hemorrhagic lesions, and with available CSF biomarkers. We visually quantified MRI-visible SVD markers (e.g., cerebral microbleeds [CMB], cortical superficial siderosis [cSS], lacunes) and their association with MRI-visible SVD features. We employed a Gaussian Mixture Model (GMM) to identify a data-driven threshold for amyloid positivity (A+). Then, we compared the prevalence of MRI-visible manifestations of SVD between subgroups applying different frameworks, namely the current MRI-based classification (probable CAA vs. mixed SVD) and a CSF biomarker-integrated classification (A+ vs. A-). Results: We enrolled 121 patients (age: 72 [66-77] years; 60% probable CAA, 40% mixed SVD with suspected CAA). The CSF A{beta}42/40 ratio showed a bimodal distribution and consistent associations with all CAA-specific radiological features. The CSF biomarker-integrated reclassification, particularly using the GMM cut-off, significantly improved the distinction between subgroups regarding CAA- and arteriosclerosis-related MRI features (e.g., cSS presence: probable CAA vs. mixed SVD: aOR=2.84 [95%CI 1.27-6.39], p=0.011; A+ vs. A-: aOR=12.68 [95%CI 4.31-37.32], p<0.001; deep lacunes presence: probable CAA vs. mixed SVD: aOR=0.20 [95%CI 0.08-0.50], p<0.001; A+ vs. A-: aOR=0.04 [95%CI 0.01-0.11], p<0.001). Notably, patients classified as A+ never demonstrated more than four deep CMBs. Discussion: A CSF biomarker-integrated classification may improve the classification of CAA compared with the current MRI-based framework. These findings are cohort-specific and would benefit from further validation, especially with a neuropathological reference. Still, these results support a future transition toward an integrated biological-radiological framework, which may refine in vivo CAA diagnosis, particularly in mixed SVD.

Investigating the conformational dynamics of intrinsically disordered proteins (IDPs) is essential to understanding how their structural heterogeneity underlies function and how their dysregulation contributes to diseases. Here, we utilized an MspA nanopore-based approach for studying the conformational dynamics and interactions of IDPs at the single-molecule level. The platform was demonstrated using the intrinsically disordered transactivation domain of tumor suppressor p53 (p53-TAD), one of the important proteins in cancer biology. We showed that MspA can stably capture p53-TAD and resolve up to six distinct current states with frequent interconversions, revealing a rich conformational landscape. The nanopore also detected the effect of a cancer-associated double mutational variant, N29K/N30D. Combining experiments with steered molecular dynamics simulations, we showed that the mutant sampled compact conformational states more frequently than wild type, consistent with previous NMR studies. Importantly, the MspA platform enabled direct monitoring of E3 ligase MDM2 binding to p53-TAD and resolved how this interaction is inhibited by anti-cancer compound epigallocatechin gallate (EGCG). Notably, EGCG stabilizes one of the six states sampled by p53-TAD, providing a mechanistic explanation for its inhibitory effect. Together, these findings demonstrate the promise of the nanopore platform for label-free monitoring of IDP conformational dynamics, modulation, binding and inhibition at single-molecule resolution.

The study of evolution between species (phylogenetics) and the study of evolution within a species (population genetics) are highly related, as the same biological mechanisms are fundamental to both fields. Although both have been studied for a long time, their joint study in a unified setting has been prevented by the different time scales they consider and the different data types they employ. A similar discrepancy holds for their whole-genome specializations, comparative genomics and pangenomics. Two active areas in these fields are genome rearrangement studies and graphical pangenomics, respectively. Since the emergence of graphical pangenomics, these have existed as separate fields, despite observations that central data structures representing genomic variants in both fields are highly similar. While there exists a wealth of theoretical results for various rearrangement models in comparative genomics, the application to pangenomic data is hampered by the limitations of rearrangement problem formulations. On the practical side, pangenomes typically contain too many individual genomes for classical problems, such as the often NP-hard parsimony problems, to be solved, or for all-vs-all comparisons using rearrangement distances to be performed. On the theoretical side, some assumptions in the formulation of rearrangement problems, such as the assumption of an underlying tree, are inadequate for many pangenomes. In this work, we propose the Complete Ancestral Reconstruction for Pangenomes (CARP) problem, which overcomes these limitations while retaining intuitive relationships to both classical rearrangement problems and pangenome graphs.

Background Understanding the mechanisms underlying brain aging and age-related pathological changes is essential for advancing brain health research. Our group previously developed a mechanistic mathematical model of healthy brain, Chamberland et al. (2024) that integrates key biological processes involved in normal aging, from which Alzheimer's disease (AD) related changes may emerge naturally. Objectives To characterize and validate this brain model by evaluating its sensitivity, calibrating its parameters, and assessing generalizability in independent populations. Methods The model represents the evolution of key biological processes associated with brain aging, including amyloid beta (A{beta}), tau pathologies, neuroinflammation, and neuronal death. After identifying the 30 most influential parameters, we calibrated the model using cognitively normal (CN) participants from the AD Neuroimaging Initiative (ADNI) database (n = 211) by minimizing a loss function composed of three outcomes (AB) plaques, tau tangles, and neuronal density). The calibrated model was then applied to the UK Biobank cohort (n = 35,899) of normal controls (aged 44-82 years). The effects of sex and APOE were evaluated using stratified simulations. Results Parameter calibration significantly reduced the prediction errors for A{beta} and tau. Neuronal density predictions showed strong agreement in the UK Biobank cohort. The variance decomposition identified APOE status as a major contributor to variability in A{beta}. Conclusion Our validated brain health model links mechanistic pathways with population data and reproduces neuronal density patterns in an independent cohort. These findings support its use as a framework for studying brain aging and investigating how Alzheimer's disease related pathological changes may emerge with aging.

Background: Serial glial fibrillary acidic protein (GFAP) trajectories have become an important framework for contextualizing evolving secondary-injury pathophysiology after moderate-to-severe traumatic brain injury (msTBI). However, total GFAP pools release and clearance signals that may be less useful for longitudinal bedside decisions than a proteoform-resolved assay. We compared total GFAP with neoGFAP, defined here as calpain-generated GFAP proteoforms intended to index active astroglial proteolysis during the subacute phase. Methods: We analyzed 651 serial serum samples from 95 msTBI patients from a previously described single-site cohort. Total GFAP and neoGFAP were measured on the same MSD platform from 6 to 240 hours after injury. Early (6 to 72 h) and late (96 to 240 h) windows, data-derived tertiles, and serial trajectory summaries were calculated directly from serial samples. Models were benchmarked against age plus admission post-resuscitation Glasgow Coma Scale (GCS) and the admission IMPACT extended risk score using five-fold stratified cross-validation. Outcomes were unfavorable outcome (GOSE 1 to 4), less-than-good recovery (GOSE 1 to 6), Disability Rating Scale (DRS) [&ge;]15, mortality, and neuroimaging worsening at 6 months. Results: The cohort contributed 95 serial biomarker profiles, with 90 participants evaluable for 6-month GOSE and 89 for DRS. Unfavorable outcome occurred in 57/90 (63.3%), and less-than-good recovery in 79/90 (87.8%). For unfavorable outcome, IMPACT plus early neoGFAP reached AUROC 0.85 versus 0.84 for IMPACT plus early total GFAP and 0.81 for IMPACT alone. For less-than-good recovery, IMPACT plus late neoGFAP achieved AUROC 0.90 versus 0.84 for late total GFAP and 0.82 for IMPACT alone. Secondary analyses for DRS, mortality, and neuroimaging worsening showed smaller differences. Conclusions: In this retrospective analysis, neoGFAP provided clearer incremental value than total GFAP for recovery-oriented monitoring, especially when late-window reassessment of patients who remained at risk for less-than-good recovery was required. Results support prospective testing of neoGFAP as a pathophysiology-informed adjunct to serial bedside decision making, repeat-assessment thresholds, and recovery stratification.

Background: On-call clinicians frequently report the anecdotal impression of 'theme shifts' during which specific acute neurological diagnoses appear to cluster. Whether such clustering reflects a statistically true and reproducible phenomenon has not been systematically investigated; the present paper examines seasonality and temporal clustering within six different acute neurological conditions. Methods: In this retrospective, single-center cohort study, we identified all patients admitted to a tertiary neurological department between July 2016 and June 2026 with acute unilateral vestibulopathy, cerebral artery dissection, generalized epileptic seizures, primary intracerebral hemorrhage, peripheral facial nerve palsy, or transient global amnesia (TGA) (n = 2,140). Monthly and seasonal distributions were assessed using chi-squared goodness-of-fit and cosinor analysis. Short-term temporal clustering was tested by Monte Carlo permutation across time windows from 24 hours to 90 days, and endogenous cluster dynamics were characterized using Hawkes self-exciting point process modeling. Results: Admissions for generalized epileptic seizures showed a statistically significant deviation from a uniform monthly distribution with a winter distribution (p<0.001 and q = 0.002), and a significant temporal clustering across time windows from 72 hours to 90 days (all q < 0.05). Peripheral facial nerve palsy presented significant clustering at the 90-day window (q = 0.029) and TGA at 60-day time window (q = 0.041) without seasonality; the diagnostic groups of acute unilateral vestibulopathy, cerebral artery dissection and primary intracerebral hemorrhage showed neither seasonality nor clustering after correction for multiple comparison. No diagnostic group showed clustering within a 24-hour window, statistically significant self-excitation in Hawkes process modelling, or a significant linear trend in monthly case counts over the study period. Conclusion: The anecdotal impression of diagnostic 'theme shifts' among on-call neurologists appears to have a measurable basis, although clustering is confined to specific conditions and rather on a time scale of weeks to months. Generalized epileptic seizures were the only diagnostic group that uniquely combined seasonality with temporal clustering, suggesting a shared trigger, while facial palsy and TGA showed episodic, yet non-seasonal clustering.